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ABSTRACT 


A  method  is  presented  which  reduces  the  number  of 
components  needed  in  a  Fialkow-Gerst  multiport  RC  trans¬ 
fer  function  synthesis.  A  relationship  is  determined 
between  the  number  of  non-zero  numerator  coefficients  in 
the  transfer  function  vector  and  the  number  of  components 
used  in  the  synthesis  of  that  vector. 
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INTRODUCTION 


Although  the  Fialkow-Gerst  synthesis  technique 

1-4 

has  been  known  for  some  time  ,  and  although  it 
has  been  extended  to  multiport  networks  by  Zeren 
and  others  '  ,  the  large  number  of  components  it  uses 

O 

is  still  a  major  problem.  Kodali  has  worked  on  re¬ 
ducing  this  number  and  an  extension  of  his  work  will 
be  presented  which  reduces  the  number  of  components 
in  three  ways.  First,  it  eliminates  components  from 
the  termination  of  the  network.  Second,  it  permits 
the  calculation  of  an  arbitrary  constant,  X,  arising 
in  the  Fialkow-Gerst  synthesis.  Finally,  the  method 
yields  transfer  functions  which  are  frequently  amen¬ 
able  to  special-case  synthesis  techniques. 

The  importance  of  RC  transfer  function  synthesis 
has  been  increased  by  the  recent  work  of  Hazony  and 

7 

Joseph  ,  which  permits  the  synthesis  of  any  RLC 
transfer  function  vector  with  one  unity  gain  ampli¬ 
fier  and  an  RC  network. 

A  computer  program  is  provided  in  the  Appendix 
which  incorporates  this  variation  of  the  Fialkow- 
Gerst  synthesis. 


CHAPTER  I 


THE  FIALKOW-GERST  TRANSFER  FUNCTION  SYNTHESIS 


1 .  Synthesis  of  a  Transfer  Function  Matrix 

The  Fialkow-Gerst  synthesis  of  a  grounded  mult 
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port  RC  transfer  function  network  is  well  known 
and  will  be  outlined  in  this  chapter  only  to  intro¬ 
duce  notation  and  provide  a  reference  for  later  dis 
cussion. 


Fig.  1.1  A  multiport  RC  network. 

We  can  describe  the  network  of  Figure  1,1  by 
the  equation 


r^outi  =  •  r^ini 


where  the  elements,  or  entries,  of 
given  by 


the  matrix  T  are 
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3- 


(1.1) 


t  .(s)  = 

qD 


a  .s^  +  a.  + 

oqi _ Ifll 


+  a 


+  b,  + 

Oq  Iq 


J3J. 


+  b 


rq 


A  .(s) 
-SI  — 

B  (s) 

q 


y  . 
-jai 

Y 

qq 


for  q  *  m+l  ,  ,  ra+n,  and  for  j  *  1  ,  ...  #  m. 

If  a  row  of  [T]  is  considered  as  a  vector  ^  then 

with  appropriate  surplus  factors,  ^  can  be  made  to 
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have  only  positive  coefficients  *  .  Each  entry  t  . 

QD 

must  then  satisfy: 


i) 


(1.2a)  ii) 
(1.2b)  iii) 


The  poles  are  distinct  and  lie  on  the 

negative  real  axis; 

0  <  2a .  .  <  b .  ^  for  q  ■  m+1,  ...  ,  m+n 

iqj  iq  ' 


D 

0  <  a.  . 
iqD 


1*0/  ...  /  r, 
j  ■  1/  ...  /  ra. 


Such  t  .  are  termed  RC  R-functions  and  a  matrix 

qD 

with  all  RC  R-function  entries  is  called  an  RC 
t-matrix. 

For  convenience  we  will  assume  throughout  this 


b  b 
o  r 


^  0. 


i 


paper  that 
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Following  Fialkow,  et  al.,  we  synthesize  only 

til 

one  row  of  an  RC  t-matrix  at  a  time.  If  that  q  row 
is  given  by  the  vector 


^^ql'  ^2' 


qrn 


of  degree*  r  with  m  vector  entries,  then  the  first 
step  in  the  synthesis  is  to  choose  a  polynomial 
(of  degree  (r-1)  with  negative  real  simple  zeros) 
such  that 


s^  +  b,  s^-^  + 

oq  iq 


+  b 


d^  ^  +  . . .  +  d  , 

Oq  r-l,q 


and 


for  all  j, 


are  RC  admittances. 

The  second  step  is  to  split  the  short  circuit 

driving-point  admittance,  Y  .  into  two  such  functions, 

qq 

Y'  and  Y"  .  Thus  for  the  first  cycle 

q<3  qq 


In  this  paper  we  are  only  concerned  with  RC 
transfer  function  vectors  whose  entries  have  a  common 
denominator.  Thus  the  term  "degree"  is  used  only  in 
reference  to  that  denominator  and  not  the  transfer 
function  matrix  as  a  whole. ^ 


1 


1 


R.  + 


^qq 


Obviously  any  method  which  accomplishes  this 
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split  is  valid.  A  method  by  Hazony  '  ,  which  will 
be  used  later,  utilizes  a  split  factor  X  so  that 


(1.3) 


■qq 


Joa!  +x(^2a!^ 

oq  ^  +  ...  +  a 


r-l,q 


+ 


+ 

q 


(i-x)  ( 


r— 1 

9a  ®  +...+Q  •  S 

Oq  ^r-2,q 

^rVr®^  +  ...  +  <3  . 

Oq  r-l,q 


where  0  <  X  <  1  but  otherwise  X  is  arbitrary  and  in 
general  it  will  have  a  different  value  for  each  syn 
thesis  cycle. 

Using  equation  (1.3)  we  can  write  the  new  ad¬ 
mittance  as 

qq 


(1.4a) 


r(i).  Oq  !  ■  _ 


**la  » 


+  b 


oq 


s 


+  d 


TTT 

r-2,q 


(1) 

r-i/q 


where 


i 
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r 


(1.4b) 


^Oq  , 


I 


(1.4c) 


b(i)  .  b.  X  +  ’  q 


iq  iq 


Oq 


r-l,q 


for  i  =  1,  ...  ,  r-1,  and 


(1.4d) 


iq 


=  d 


iq 


bf^^d  . 

■  >q-..f-i./q 
J(i) 

^r-l,q 


for  i  *  0,  ...  ,  r-2. 


(2) 


Similarly  we  can  write  Y  as 

qq 


(1.5a) 


r(2) 

qq 


,  (2)  r-l  (2)  r-2  (2) 


Iq 


r-l,q 


where,  for  i  =  1,  ...  ,  r-l, 


(1.5b) 


b,!!'  =  b,„(i-x)  -  iq - +  ■^q  ^-^.q  , 

Oq  “^r-l.q 


iq  iq 


(1.5c) 


b^^^  =  b  ,  and 
rq  rq  ' 


(1.5d) 


d!^^  =  d. 
iq  iq 
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The  resistor  and  capacitor  removed  are  given  by 


(1.6a) 


C,  =  /  d  , 

1  r-1  '  r-1, 


farads  and 


(1.6b) 


d  /  b|^^ 
oq  '  Iq 


ohms, 


These  equations  emphasize  the  adaptability  of 
this  synthesis  method  to  a  computer  solution. 

The  final  step  in  the  synthesis  cycle  is  the 
computation  of  the  two  reduced  transfer  function 
vectors  such  that  each  vector  entry,  as  well  as  the 
sum  of  the  entries,  is  an  RC  R-function.  Hence,  for 
the  first  cycle  _t  becomes 


-q 


B 


A<1) 

am 


and 


A<?> 


(2) 


qm 


B 


where  the  numerator  coefficients  satisfy 


(1.7a) 


,  a  .  .  a'^) 

Oq:  Oq:  rq:  rq: 
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(1.7b) 


a,  .  = 

iq3  iq3  iq3 


for  i  =  1,  ,  r-1,  and  all  j, 

and  the  denominator  coefficients  satisy 


(1.7c) 


b  *  b(l)  b  = 

Oq  Oq  '  rq  rq 


(1.7d) 


,  _  .  (1)  .  ,  (2)  .  1 

‘’iq  -  ^  ‘’iq  '  ‘ 


..  ,  r-1. 


Equation  (1.7)  can  be  represented  as 


*  s  A 


(1) 

■qj 


+  A 


(2) 

■qj 


for  all  j,  and 


B 


»  s  B 


(1) 


+  B 


(2) 


Obviously,  to  preserve  the  RC  R-f unction  charac¬ 
teristics,  the  new  coefficients  must  satisfy  equation 
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(1.2).  Hazony  '  has  introduced  a  proportional 
method  for  calculating  the  reduced  transfer  function 
vector  numerators  which  uses  the  following  two 
equations : 


(1.8a) 


t 


for  i  =  1,  ...  ,  r-1,  and  all  j. 
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and,  by  equation  (1.7a); 


(1.8b) 


a<2»  .  a  . 

rq:  rqj 


for  all  j. 


Another  method  of  calculating  the  transfer  functions 
will  be  presented  in  Chapter  2. 

After  one  cycle  the  vector  network  looks  like 
Figure  1.2;  where  is  a  network  with  a  transfer 

function  vector  s  driving-point  admittance 

yd).  p(2)  g  similar  network. 

qq  q 


Pig.  1,2  The  results  of  one  synthesis  cycle. 

The  above  process  is  repeated  until  unity  degree 


transfer  functions  are  obtained;  at  which  time  the 
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degree  three  with  two  vector  entries. 
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functions  are  synthesized  directly  and  the  network 
will  appear  as  in  Figure  1.3.  In  the  subsequent  dis¬ 
cussion  we  will  refer  to  the  set  of  resistors  and  cap¬ 
acitors  obtained  from  splitting  the  driving-point 
admittance  as  the  tree  of  the  network  and  the  com¬ 
ponents  obtained  from  the  unity  degree  transfer 
functions  as  the  network  termination. 

Once  all  of  the  row  vectors  have  been  synthesized, 

all  of  the  inputs  are  connected  in  parallel  to 

til 

yield  one  common  j  input  to  the  whole  network.  As 
there  is  only  one  output  there  is  no  interconnec¬ 

tion  between  the  output  terminals.  The  interconnec¬ 
tion  takes  the  form  of  Figure  1.4. 


1  o- 


D  O- 


m 


m 


m+1 


m+1 


■^+n 


m+n 


m+n 


Fig.  1.4  Interconnection  of  synthesized  vectors. 
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2 .  The  Auqumented  Transfer  Function  Vector 

When  synthesizing  the  unity  degree  transfer 
function  vectors  the  grounded  components  are  values 
proportional  to  the  differences  between  the  denomin¬ 
ator  coefficients  and  the  sums  of  the  numerator  co¬ 
efficients.  It  is  frequently  convenient  to  make  this 
"transfer  function"  to  ground  explicit:  the  comple¬ 
mentary  transfer  function^.  t^Q,  is  given  by  t^^  = 

1  -  S  t  .  and  the  augmented  transfer  function  vector, 

t,  is  the  vector  t  with  the  entry  t^^  added,  i.e., 

QO 

^ql'  ***  '  ^qm^* 

Since  the  augmented  vector  has  the  property  that 


(1.9) 


m 


=  1  , 


there  are  no  components  to  "ground"  when  is  synthe- 
sized.  Instead,  there  are  components  to  the  0^  input 
port  and  that  port  has  a  voltage  generator  of  zero 
volts  connected  to  it. 

Let  us  briefly  examine  the  effect  of  inter¬ 
changing  one  of  the  input  terminals  with  the  connec¬ 
tion  to  ground.  Using  Figure  1.1,  Figure  1.5a  can 
be  described  by  the  equation 


Sl"lO  + 


+  t  .E.,,  + 
qi  dO 


+  t  E  . 
qm  mO 
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Subtracting  EjQ  from  both  sides  and  then  adding  and 

subtracting  t^^E^Q  on  the  right  side  yields  (recall 

that  E .  .  *  E .  .  -  E . ^) 

13  lO  3O 


E  . 

qi 


t  ,E,  .  +  t  -E. .  + 
ql  I3  q2  2d 


+  (t  .-1-t  .)E.^ 

q:  qD  30 
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which  describes  the  transfer  function  vector  for  the 

th 

new  network  of  Figure  1.5b  in  which  the  j  terminal 
is  common. 

It  is  apparent  that  an  exchange  of  the  j  input 

terminal  with  ground  amounts  to  interchanging  t^Q 

with  t  .  in  t"*",  which  is  really  no  more  than  renaming 
q3  -q 

the  terminals.  This  will  be  used  later  to  justify 
arbitrary  rearrangements  of  the  augmented  tr^insfer 
function  vector  entries. 

When  t  .  is  zero,  there  is  no  connection  to  the 

qi 

j  terminal  and  therefore  the  vector  can  be  regarded 
as  an  (m-1)  entry  vector.  This  fact  will  be  utilized 
in  Section  4  of  Chapter  2. 


CHAPTER  2 


A  METHOD  OF  SPLITTING  AN  RC  TRANSFER  FUNCTION  VECTOR 
BY  CREATING  DEGENERACIES 

1.  Introduction 

0 

Kodali  has  developed  a  method  of  calculating  the 
reduced  transfer  functions  t*  and  t"  which  utilizes 
both  the  arbitrary  splitting  factor  \  of  equation 
(1.3)  and  the  arbitrariness  of  the  polynomial  D^. 

The  method  outlined  below  extends  his  results  to 
vector  transfer  functions  and  explicitly  considers 
functions  of  arbitrary  degree.  The  method  also  pro¬ 
vides  useful  values  for  the  split  factor  "K. 

As  we  saw  in  the  last  chapter,  the  actual  method 
used  to  obtain  the  vectors  _t‘  and  is  arbitrary 
provided  that  the  coefficient  conditions  are  satis¬ 
fied.  The  method  presented  below  forces  the  numer¬ 
ator  coefficients  to  satisfy  the  equalities  of  equa¬ 
tion  (1.2)  as  much  as  possible.  This  leads  to  fewer 
components  being  needed  in  the  synthesis  of  the  unity 
degree  transfer  functions. 

2.  Transfer  Function  Calculations 

The  method  below  considers  the  coefficients  of 

r— i 

s  as  a  group  of  numbers  which  can  be  split  into 
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two  new  groups  such  that  the  sum  of  one  of  the  new 


groups  is  either  zero,  b*  ,  or  b".  One  of  the  new 
groups  is  then  assigned  to  _t'  and  the  other  to  as 
the  numerator  coefficients  of  s^  ^  ^  and  s^  ^  re¬ 
spectively.  The  reader  can  best  follow  the  general 

method  of  calculating  t'  and  f  by  considering  a 

-q  -^3 

particular  case  first. 

Example  2.1,  Suppose  that  the  numerator  coefficients 

3 

of  s  in  a  five  entry  transfer  function  vector  of  de¬ 
gree  seven  are  given  by  =  3,  a^2  “  5,  a^^  “  2, 
a^^  *  9,  and  a^^  =  1,  or  more  compactly,  as 
^  =  (3,  5,  2,  9,  1),  and  chat  the  denominator  co¬ 
efficient  is  b^  =  30.  Suppose  also  that  in  splitting 
the  admittance  by  using  equations  (1.4)  and  (1.5)  we 
obtain  b^  =  21  and  b^  *  9.  Then  using  the  proportion¬ 
al  method  of  equation  (1.8)  for  calculating  the  reduced 
transfer  functions  we  obtain 


7,  7,  63,  _J_. 
2  5  10  lO' 


and  a ' 


3,  3,  27,  _3v 
UO  2  5  10  10^ 


In  splitting  these  coefficients  we  only  require 
that  equation  (1.2)  be  valid,  namely  that 


i 
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and  each  coefficient  is 


a .  +  a "  .  =  a 


It  is 


4j  “4j 
are  also  met  by  each  of 


non-negative,  and  that 
apparent  that  these  conditions 
the  following: 


i) 

^  =  (3, 

5, 

2, 

0, 

II 

O 

0, 

0, 

9, 

0), 

iia) 

o 

II 

0, 

1, 

9, 

1),  =  (3, 

5, 

1, 

0, 

0), 

iib) 

o 

II 

0, 

2, 

9, 

0),  =  (3, 

5, 

0, 

0, 

1), 

iii) 

=  (3, 

5, 

2, 

9, 

II 

0, 

0, 

0, 

0)  . 

The  advantage  of  this  second  type  of  calculation 
lies  in  the  number  of  zero  coefficients  introduced, 
since,  in  the  last  stage  of  a  synthesis,  every  non¬ 
zero  coefficient  in  the  augmented  transfer  function 
vector  is  proportional  to  the  inverse  value  of  a 
resistor  or  capacitor.  In  the  following  pages  it  will 
be  shown  that  one  can  always  calculate  ^  and  ^  in 
this  manner  and  that  the  worst  possible  case,  in 
terms  of  the  number  of  non-zero  coefficients,  yields 
four  instead  of  five  zeros  for  a  five-entry  vector. 

For  the  sake  of  clarity  let  us  drop  the  sub- 

scripts  q  and  i,  i.e.,  let  ^  ^  and  let  = 

a.s^-\ 

D 

The  three  basic  splits  of  Example  2.1  correspond 
to  the  three  methods  that  can  be  used  to  make  the  co¬ 
efficients  equal  to  0,  b'  ,  or  b”.  Assuming  that 
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I 

I 

b'  <  b",  the  coefficients  of  s^”  in  the  numerator 

and  denominator  of  the  vector  ^  must  conform  to  at  ^ 

least  one  of  the  following  three  classifications: 

Class  i)  a^  >  b'  for  some  k  between  1  and  m. 

With  a^  >  b*  we  can  split  a^^  into 
two  parts;  the  first  will  equal  b'  ;  the  re¬ 
mainder  and  all  the  other  coefficients  will 
form  _t " . 

Putting  this  split  in  equation 
form  we  have 

(2.1a)  a'^  =  b'  and  a"^  =  a^-  b' , 

(2.1b)  a'  .  =  0  and  a".  =  a.  ,  for  all  j  k. 

D  i  j 

It  follows  that  such  a  split 
satisfies  all  of  the  coefficient  conditions 
and  that  it  creates  (m-1)  zero  coefficients. 

If  a^  >  b"  we  may  wish  to  use  b" 
and  thus  rewrite  equation  (2.1)  with  '  and 
"  interchanged. 

Class  ii)  a^  <  b'  for  every  j  and  yet 
S  a .  >  b". 

j  ^ 

Since  no  one  coefficient  is  large 
enough,  a  partial  sum  must  be  formed  such 

I 

that  by  adding  a  fraction  of  a^  to  it,  the  I 


I 
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partial  sum  is  equal  to  b”.  The  re¬ 
mainder  of  aj^and  the  unused  coefficients 
will  form  _t'  .  Thus  we  pick  k  such  that 

k-1  k 

(2.2)  2  ^  a.>b" 

j=i  3  j=i  3 

and  calculate  the  new  coefficients  by 
setting 

(2.3a)  a^  =  0  and  a''^  =  a^  for  j  =  1/  ...  /  k-1, 

k  k-1 

(2.3b)  a'  =  -  b”  +  2  a .  and  a"  =  b"  -  S  a.  , 

j=l  3  K  j=l  3 

(2.3c)  a^  =  a^  and  aV  =  0  for  j  *  k+1,  ...  ,  m. 

Since  the  entries  of  a  transfer 
function  vector  can  be  moved  about  in  the 
vector,  one  is  free  to  group  the  coefficients 
in  any  manner.  At  least  (m-1)  coefficients 
will  be  zero  but  sometimes  a  careful  group¬ 
ing  will  lead  to  a  sum  which  equals  b”, 
making  a^  zero  (cf.  Example  2.1  iib)  . 

Since  b'  <  b",  one  could  write  a 
set  of  equations  which  assign  a^  with  re¬ 
spect  to  b*  .  These  equations  would  merely  be 
equations  (2.2)  and  (2.3)  with  '  and  ”  inter¬ 
changed  . 
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It  is  easily  shown  that  the  co¬ 
efficient  conditions  are  satisfied  by 
equation  (2.3):  the  aV  were  constructed  to 
satisfy  the  conditions  with  respect  to  b"; 
summing  the  a^  we  obtain 

m  km 

7.  a'.  »  (0)  +  (-b"  +  I  a.)  +  (  la.) 
j=l  ^  ^  j=k+l  ^ 


m 

=  ^  a .  -  b” 

j=i  ^ 


but  since  I  a.  ^  b  and  since  b  -  b”  =  b'  it 

j  ^ 

is  apparent  that  the  a^  do  indeed  satisfy 
the  coefficient  conditions  with  respect  to 
b'  . 


Class  iii) 


I  a.  <  b". 

j  ^ 

Under  this  condition  there  is  no 


need  to  split  coefficients  since  obviously 
the  a^  already  satisfy  the  coefficient  con¬ 
ditions  with  respect  to  b"  and  thus  can  be 
put  directly  into  _t",  i.e., 


(2.4) 


a'  .  *  0  and  a'.'  *  a.  for  all  j. 
3  3  3 


This  creates  m  zero  coefficients. 
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Of  course,  if  the  sura  is  less  than 
or  equal  to  b'  we  are  free  to  use  equation 
(2.4)  with  '  and  "  interchanged. 

Regardless  of  which  of  the  above  three  methods 

X  0 

is  used,  the  coefficients  of  s  and  s  must  be  treated 
as  in  equation  (1.7a). 

When  the  transfer  function  is  a  scalar,  the  more 
convenient  classifications  are:  1)  a^^  <  b”  to  use 

method  iii)  and  2)  a^>  b”  to  use  method  i)  .  The 

second  vector  classification  for  method  ii)  vanishes 
trivially. 

3 .  The  Number  of  Components  in  a  Synthesis 

To  evaluate  the  effect  of  the  above  split  on  the 
number  of  components  needed  in  a  synthesis  two  de¬ 
finitions  will  be  introduced. 

Definition  1.  A  numerator  coefficient,  a. .,  of  t  is 
-  -  —  — 

said  to  be  degenerate  if  it  is  zero. 

Definition  2.  The  degeneracy,  6  ^  .  of  a  transfer 
function  vector  t  is  the  number  of  degenerate  co¬ 
efficients  in  its  augmented  vector  t"*". 

In  a  unity  degree  augmented  transfer  function 
vector  each  non-zero  coefficient  is  a  component  and 
thus  each  degeneracy  is  a  missing  component  (see 
Figure  2.1).  With  this  in  mind,  we  will  proceed  to 


calculate  the  number  of  degeneracies  at  the  end  of 
a  synthesis.  Once  this  number  is  known,  determining 
the  number  of  components  is  a  trivial  matter. 


missing 


X (n)_  r  s  1  s-t-2  ~| 

Fig.  2.1  Synthesis  of  ^  ”  [2s+3  '  2s+3  '  2s+3  J 

•  ..•L  ..(n)  2s+3  j  r  (n) 

with  Y  =  — ^ —  and  6^  -2. 

In  all  of  the  discussion  and  theorems  that 
follow  it  is  assumed  that  one  method  of  calculating 
the  reduced  transfer  functions  is  used  exclusively  in 
a  synthesis. 
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Theorem  3.  The  degeneracy  of  a  transfer  function 
vector  of  degree  r  with  m  entries  satisfies  6^  <  in(r+l). 

Proof .  In  the  vector  _t^  there  are  (m+l)  (r+1) 
numerator  coefficients  of  which  at  least  (r+1)  are 
non-zero  by  equation  (1.9).  Hence  there  are  at  most 
m(r+l)  degeneracies.  To  exemplify  the  extremes: 

=  0  for  t  =  [^f^] 

6  .  =  6  for  ^  =  \-j — -  ,  -= - — - 

[s  +  3s  +  1  s^  +  3s  +  1  . 

QED 

Notice  that  the  degeneracy  of  a  network  is  the 
sum  of  the  degeneracies  of  the  vectors  which  describe 
that  network  and  as  such  it  will,  in  general,  in¬ 
crease  as  the  synthesis  proceeds. 

Theorem  4.  Using  the  degeneracy  method  of  splitting 
t,  the  net  increase  in  the  degeneracy  of  a  network  is 
between  m(r-l)  and  (m+1) (r-1)  at  each  cycle,  i.e., 

6^^^+  6^  +  d,  where  m(r-l)  <  d  <  (m+1) (r-1) . 

Proof .  Upon  examining  the  degeneracy  split  we 

see  that  only  (r-1)  coefficients  are  involved  in  the 

r  0 

calcula  tion  since  the  coefficients  of  s  and  s  are 


■j 

i 
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not  split.  Each  of  these  (r-1)  coefficients  intro¬ 
duces  either  m  or  (m+1)  new  degeneracies,  depending 
upon  whether  or  not  the  remainder  is  zero.  Hence  the 
asserted  limits. 

QED 

Theorem  5 .  Using  the  degeneracy  split,  if  6 =  A  for 
a  vector  t ,  then  the  final  network  will  have  (A+D) 
degeneracies,  where 

m(2^-  r-1)  <  D  <  (m+1) (2^-  r-1)  . 

Proof .  Using  Theorem  4,  after  the  first  cycle 
6|^^+  6^^^=  A  +  ,  where  m(r-l)  <  <  (m+1)  (r-1). 


For  the  second  cycle 


(3) 


+  6 


(4) 


+  6 


(5) 


af  =  A 


+  d,.  +  2d. 


where  m(r-2)  <  <  (m+1)  (r-2)  since  r  was  decreased 

by  one.  We  have  used  the  notation  of  Figure  1.3  that 
_t^^^  splits  into  and  _t  . 

Similar  terms  are  added  for  the  (r-1)  times  that 
the  transfer  function  is  split  to  reduce  the  degree 
nr- 


iho  final  result  is 
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A  +  in(r-l)  +  2m(r-2)  +  4in(r-3)  +  ...  +  2^ 


< 


2^-2 


(i) 

+ 


< 


A  +  (m+l  )(r— 1)  +  2  (m+l )  ( r-2 )  +  ,..+2^^  (m+l )  , 


which  is  equivalent  to 


A  +  m(2^-r-l) 


< 


(i) 

+ 


<  A  +  (m+l) (2^-r-l) . 


QED 

In  Example  2.2  below  the  results  of  the  computer 
program  in  the  Appendix  are  used  to  show  how  the 
number  of  degeneracies  increases  as  the  degree  of  the 
transfer  function  vectors  decreases. 

Example  2.2.  Suppose  that  the  transfer  function 
vector  is  given  by 


8.88s^+18s^+17.76s  s^+5s^+4 

s^+9s^+23s^+18s+4  '  s^+9s^+23s^+18s+4 


and  the  driving-point  admittance  by 


33 


s"^  +  9s^  +  23s^  +  18s  +  4 
4s^  +  27s^  +  46s  +  18 
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Notice  that  =  7  so  that  by  Theorem  5  there 
will  be  between  29  and  40  degeneracies  in  the  final 
vectors.  After  the  first  cycle  there  are  13  de¬ 
generacies  since 


s^+5s 


7.3ls^-l-9.25s+6.14  _ 

s^+7.43s'^  1 14.25S+6.14  '  s^+7 . 43s^+14 . 25s+6 . 14 


1.57s^->-8.75s^-H1.6s  _ 4 _ 

1.57s^+8.75s^+11.9s+4  '  1 . 57s^+8 . 75s^+ll . 9s+4 


The  second  cycle  (with  "K  =  0.5  everywhere  except 
for  splitting  of  »  when  it  is  0.387)  gives 

10  more  degeneracies  in  the  vectors 


^3 


(3)  _ 


5.15s 


s^+5 


s‘+5.27s+5 


s‘+5.27s+5 


^3 


(4)  _ 


2.16s  +9. 258+6.14 
2.16s^+9.25s+6.14 


-3 


(5)  _ 


1.57s  +5.80S+3.19 
L1.57s^+5.80s+3.19 


0 


^3 


(6)  _ 


2.95s  +8. 43s 


.2.95s^+8.67s+4  '  2 . 95s^+8.67s+4 J 


The  third,  and  final,  cycle  gives  the  following  unity 
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degree  transfer  function  vectors: 

,(7)  ,|-  3.51 


(8) 


(9) 


^3 


.  (10) 
-3 


.  (13) 
-3 


.  (14) 
-3 


s  +  3.62 

1.64s 

'  s  +  3.62 

,  5 

1.64s  + 

5 

1.64s  +  5 

2.16s  + 

4.93 

2.16s  + 

4.93 

1  ^ 

4.32s  + 

6.14 

A 

4.32s  + 

6.14 

1.57s  + 

3.47 

A 

1.57s  + 

3.47 

1  u 

2.33s  + 

3.19 

A 

2.33s  + 

3.19 

#  0 

2.95s  + 

2.69 

A 

2.95s  + 

2.69 

#  0 

5.74s 

4 

] 

] 


5.98S+4 


] 


The  network  which  synthesizes  given  in 

Figure  2.2  and  since  the  final  vectors  have  a  total 
of  30  degeneracies  there  are  only  32  components  in 
the  network. 

Returning  to  the  task  of  calculating  the  number 
of  components  in  a  network  we  have: 

Theorem  6.  The  maximum  number  of  components  needed 
in  a  Fialkow-Gerst  synthesis  employing  only  the  de- 


I 
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Fig.  2.2  Computer  solution  of  Example  2.2. 
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r  4*  1 

qeneracy  split  is  (2  -24-in (r+l)  )  . 

Proof .  Considering  the  absolute  worst  case  of 

jr 

D  =  0,  by  Theorem  5  there  are  at  least  m(2  -r-1)  de¬ 
generacies  when  the  synthesis  is  complete.  Since 
there  are  2  (m+1)  coefficients  among  the  unity  degree 
transfer  functions,  -he  difference,  or  2  +in(r+l),  will 
appear  as  components  in  the  network  termination.  Add- 
ing  to  this  number  the  (2  -2)  components  in  the  net¬ 
work  tree  we  obtain  the  asserted  result  of  2+m(r+l) 

components . 

A  particular  transfer  function  vector  may  use 
fewer  components  or  a  practical  realization  may  use 
more  but  ( 2*^^^-2+m  (r+1))  is  the  ''  maximum  number  of 
components  needed”,  i.e.,  it  is  a  sufficient  number. 
Corollary  7.  The  maximum  number  of  componen^ts  nefeded 
in  a  Fialkow-Gerst  synthesis  employing  only  the  pro¬ 
portional  split  is  ( 2^^^-2-fm2^)  . 

Proof.  It  is  evident  that  if  the  transfer 
function  vector  has  no  degeneracies  at  the  start 
then  the  proportional  split  will  not  introduce  any. 
Hence  none  v/ill  appear  in  the  unity  degree  transfer 
function  vectors  and  all  the  termination  components 

will  be  present.  These  termination  components  and 

r4*  1  r 

the  tree  components  add  to  (2  -2+m2  ). 

QED 
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The  major  difference  between  the  two  methods  is 
the  factor  of  (r  +  1)  instead  of  2  in  the  maximum.  The 
difference  between  the  methods  in  terms  of  the  minimum 
number  of  components  obtainable  is  not  as  dramatic. 

In  fact,  the  real  difference  lies  not  in  the  number 
of  components  but  in  the  number  of  inputs,  or  vector 
entries,  that  each  method  can  tolerate  before  the 

IT  1 

absolute  minimum  of  (2  -2)  becomes  unobtainable. 

It  will  also  be  sean  that  the  degeneracy  split  permits 
a  larger  variety  of  minimal  vector  forms. 

Lemma  8 .  The  minimum  number  of  components  obtainable 

JT  “^1 

using  the  proportional  split  for  m  <  (r+1)  is  (2  -2) 

r+l 

and  for  m  >  (r+1)  is  (2  -2)  +  (m-r-1) . 

Proof .  Regardless  of  the  method  of  calculation, 
r- 1 

2  unity  degree  transfer  function  vectors  will  be 
produced  in  a  synthesis.  Since  each  of  these  has  at 
least  two  non-zero  coefficients  in  its  augmented  vector 
there  must  be  at  least  2  components  in  the  network 
termination.  When  the  components  in  the  network  tree 

JT  1 

are  added  there  is  a  total  of  (2  -2)  components. 

Since  a  maximally  degenerate  vector  will  use  a 
minimum  number  of  components  the  problem  becomes: 
when  is  it  impossible  to  use  a  maximally  degenerate 
vector?  For  the  proportional  split  this  cutoff 
occurs  when  all  the  denominator  coefficients  have 
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been  used  to  form  entries,  i.e.,  when  m  >  (r+1) .  If 

the  extra  entries  are  created  by  using  fractions  of 

0  r 

the  coefficients  of  s  or  s  there  will  be  an  increase 
of  only  one  component  for  each  of  the  fractional  en¬ 
tries.  This  is  due  to  the  fact  that  the  coefficients 
0  r 

of  s  and  s  are  not  split  in  the  synthesis. 

The  existence  of  the  vectors  is  demonstrated  by 

2  2  " 

^ _  3s  1  s^ 

~  2s^+3s+l  *  2s^+3s+l  2s^+3s+l  2s^+3s+l 

which  only  requires  seven  components. 

It  should  be  noted  that  if  any  other  coefficient 
is  used  then,  as  the  vector  is  split,  this  extra  co¬ 
efficient  will  appear  in  more  and  more  vectors  until 

r  0 

it  becomes  the  s  or  s  coefficient.  These  excess 
coefficients  will  of  course  produce  excess  components. 
Theorem  9.  The  minimum  number  of  components  obtain- 
able  using  only  the  degeneracy  split  for  m  <  2  is 
(2^'*‘^-2)  and  for  m>2^  is  (2^'*‘^-2)  (m-2^)  . 

Proof.  The  calculation  is  the  same  as  in  Lemma  8 
and  the  transfer  function  vector  forms  are  essentially 
the  same.  The  entry  limit  of  2  comes  from  the  fact 
that  a  number  which  is  in  the  original  vector  numerator 
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can  be  carried,  unaltered,  to  one  of  the  unity  degree 
transfer  function  vectors.  Thus,  corresponding  to 
Example  2.2,  by  starting  with  the  vector 


1  r  4 

1 - 3 - 2 -  ® 

s  +9s  +23s  +18S+4 


1.57s^  ,  2.95s^,  2.33s^ 
3.47s^  ,  6.14s  ,  4I 


one  would  get  terminal  vectors  such  as 


[s+3.62  ' 

Fo  _ L.57s_ 

1.57S+3.47 


0, 


0, 


3A1 

s+3.47 


and  this  network  would  require  only  the  minimum  thirty 
components.  Obviously  the  limit  is  the  number  of  de- 
nominator  coefficients,  2  ,  that  are  available  in  the 
termination  vectors,  not  the  (r+1)  in  the  original 
vector,  as  in  the  proportional  split. 

When  m  >  2  fractions  of  coefficients  may  be  used 
to  obtain  a  minimum  of  components  but  they  are  no 
longer  restricted  to  any  particular  coefficient. 

QED 


Applying  Theorems  6  and  9  to  Example  2.2  the 
synthesis  will  yield  between  thirty  and  forty  compo¬ 
nents  for  the  degeneracy  method  and  between  thirty  and 
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sixty-two  for  the  proportional  method.  After 
synthesizing  it  both  ways,  it  is  found  that  the  pro¬ 
portional  split  synthesis  used  forty-four  components 
and  the  degeneracy  split  used  thirty- two.  In  the 
next  section  it  will  be  shown  that  even  the  latter 
number  can  be  reduced. 

4 .  Hybrid  Synthesis 

The  degeneracy  split  greatly  enhances  the 
possibility  of  obtaining  special  transfer  functions 
which  can  be  synthesized  by  non-Fialkow-Gerst  methods, 
at  a  savings  of  components.  The  example  below  illus¬ 
trates  this  point. 

Example  2.3.  Examining  the  transfer  functions  of 

(4)  (5) 

Example  2.2  we  see  that  bo  th  _t ^  and  ^t^  are  trivial 

functions  of  the  form  t  =  [1»0],  with  second  degree 
driving-point  admittances  given  by  the  computer  as 

2.16s^+9.25s+6.14 
0.65S+1.19 

1.57s^+5.80s-l-3.19 
0.86S+1.58 


and 


.(4)  ^ 
■33 


.(5)  ^ 
•33 


Such  transfer  functions  can  be  synthesized  as  in 
Figure  2.3  by  a  ladder  network.  When  this  is  done 
for  the  problem  at  hand  we  obtain  Figure  2.4,  which 
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has  only  twenty-eight  components — two  less  than  the 
absolute  minimum  that  is  obtainable  with  the  pure 
Fialkow-Gerst  synthesis. 


1  O- 


“S 

-O  3 


20 


Fig.  2.3  Synthesis  of  ^  =  [1^0] 


5 .  Optimization  of  the  Degeneracy  Split  Using  the 
Split  Factor  'h 

The  degeneracy  split  is  based  upon  the  idea  of 
splitting  a  numerator  coefficient  such  that  the  new 
coefficients  equal  0,  b'  «  or  b".  In  general,  this 
leaves  a  non-zero  remainder  for  the  other  vector. 

Let  us  now  attack  the  problem  from  the  other  direction 
by  forcing  the  denominator  coefficients  to  be  equal 
to  the  numerator  coefficients. 

Recall  that  the  new  vector  denominators  are 
calculated  in  equation  (1.3)  when  the  admittance  is 
split.  Recall  also,  that  in  splitting  the  admittance 
we  use  a  split  factor  >  which  is  arbitrary  but  re¬ 
stricted  to  the  range  of  numbers  between  zero  and 


I 
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Fig.  2.4  Hybrid  synthesis  of  Example  2.2 
(modification  of  Fig.  2.2). 
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one .  Since  equations  (1.4c)  and  (1.5b)  relate  to 

"K  they  can  be  used  to  calculate  X  for  given  values  of 

b.  . 

1 

Several  values  of  b^  are  useful  in  that  they  in¬ 
crease  the  degeneracy  of  the  transfer  function  vector. 
For  example,  in  equation  (2.1)  we  would  like  the  re¬ 
mainder  a"^j^  to  be  zero,  hence  i  ~ 

If  a"^^  =  a|^^=  0  then  equation  (1.4c)  will  yield 


(2.5a) 


X  = 


a  -  b-d  ./d_ 
ik  0  i'  0 


b.  -  b^d^/dQ  -  b^d._^/d^_^ 


( 2 ) 

If  a'V^  =  a^^  =  0  then  equation  (1.5b)  gives 


(2.5b) 


X  = 


a.,  +  b-d./d-  -b. 

_ ik  0  r  0  1 

b.  -  hfy.d./df.  -  b  d.  ./d 

1  0  i'  0  r  i-r  r-1 


Similarly  in  equation  (2.4)  we  would  like  a' ^ 
to  be  zero  and  this  can  be  done  by  using  equation  (2.5) 


with  y,  a.  .  substituted  for  a 

j 

also 


.  Equation  (2.3)  can 
k-1 


be  used  by  substituting 


equation  (2.5) . 


j=l 


a .  .  for  a in 
13  ik 


*  8 

Kodali  has  investigated  the  conditions  that  an 
RC  transfer  function  must  meet  if  X  equals  zero  or 
one . 
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6  .  Sununarv  and  Conclusion 

In  this  chapter  we  have  shown  the  advantage  of 
the  degeneracy  split  over  the  proportional  split,. 
Example  2.3  showed  one  of  these,  namely  the  creation 
of  easily  realized  transfer  functions.  This  use  of 
the  degeneracy  split  appears  to  be  promising  and  will 
be  the  subject  of  future  investigations. 

Secondly,  the  use  of  a  calculated  value  of  the 
split  factor  Tv  is  important  since  the  introduction 
of  even  one  degeneracy  early  in  the  synthesis  of  a 
vector  leads  to  many  additional  degeneracies  in  the 
unity  degree  transfer  functions. 

Finally,  the  most  significant  contribution  of 
the  degeneracy  split  is  that  it  can  be  used  on  any 
RC  transfer  function  vector,  of  any  degree,  and  that 
it  greatly  reduces  the  number  of  components  needed  to 
synthesize  that  vector.  Table  2.1  shows  how  large 
this  difference  can  become  for  transfer  functions  of 
low  degree.  The  fact  that  these  maxima  were  predicted 
by  counting  a  particular  type  of  coefficient  suggests 
that  other  properties  of  networks  may  also  be  pre¬ 
dicted  from  the  original  coefficients. 


m\ 

2 

3 

4 

5 

10 

1 

9/10 

18/22 

35/46 

68/94 

2057/3070 

2 

12/14 

22/30 

40/62 

74/126 

2068/4094 

■5 

15/18 

26/38 

45/78 

80/158 

2079/5118 

4 

18/22 

30/46 

50/94 

86/190 

2090/6142 

5 

21/26 

34/54 

55/100 

92/222 

2101/7166 

Table  2.1 
components 
or  the 


A  comparison  of  the  maximum  number  of 

in  a  synthesis  using  the  degeneracy  split, 

proportional  split,  M  . 

P 


Appendix 


COMPUTER  SYNTHESIS 

A  computer  program  has  been  written  for  syn¬ 
thesizing  an  RC  transfer  function  vector  and  its 
driving-point  admittance  using  the  Fialkow-Gerst 
method.  The  program  incorporates  several  of  the 
techniques  introduced  in  Chapter  2  and  allows 
for  hybrid  synthesis  by  printing  all  the  calculated 
vector  and  admittance  coefficients.  For  networks  with 
only  one  or  two  inputs,  a  subprogram  is  available 
for  scaling  the  components  to  a  given  frequency  and 
for  placing  the  component  values  in  their  appropriate 
place  in  a  printed  circuit  board. 

A  general  flow  chart  of  the  program  is  given  in 
Figure  A.l  and  it  is  assumed  that  the  reader  is  suf¬ 
ficiently  familiar  with  ALGOL  60  that  the  details  will 
be  evident  from  the  program  itself.  The  one  notation- 
al  deviation  from  the  text  arises  from  a  programming 

problem - when  the  admittances  and  transfer  function 

vectors  are  split,  the  even  numbered  parts  have  sub¬ 
scripts  which  differ  from  the  text  by  unity,  thus 
A(I,P,J)  =  A(I,  2P+1,  J)+A(I-1,  2P+2,  J)  instead  of  a^^  * 
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Print 
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Fig.  A,1  General  flow  chart 
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The  great  flexibility  of  most  synthesis  procedures 
creates  problems  when  an  attempt  is  made  to  write  a 
program  for  them  and  the  ELalXow-Gerst  is  no  exception. 
Some  of  the  design  decisions  are  mentioned  below  so 
that  the  reader  will  be  aware  of  the  program's  short¬ 
comings  . 

When  several  split  factors  are  available  for  a 
vector  the  program  picks  the  one  nearest  to  0.5  in 
order  to  minimize  the  range  of  component  values.  The 
program  makes  no  attempt  to  consider  future  vectors* 

Using  the  chosen  X,  the  admittance  is  split  and 
the  resistor  and  capacitor  are  removed.  When  the 
transfer  function  vectors  and  ^^2p+2) 

calculated  a  set  sequence  is  followed.  If  method  i) 
is  used  then  the  largest  numerator  coefficient,  MXA, 
is  used  as  a^  and  the  smallest  denominator,  MNB,  is 
used  as  b'  in  equation  (2.1).  This  was  done  to  mini¬ 
mize  the  range  of  component  values.  If  method  ii) 

is  used,  no  special  grouping  or  testing  is  done - 

the  coefficients  are  summed  in  order,  starting  with 
entry  1. 

For  these  and  other  reasons  the  program  is  sub¬ 
optimum  and  it  is  suggested  that  the  user  try  all 
possible  permutations  of  the  vector  entries  in  the 
input  data,  including  using  the  complementary  term 
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as  an  input.  Using  the  complementary  transfer 
function  as  an  input  will  also  emphasize  the  round¬ 
ing  errors  of  the  computer  which  appear  as  grounded 

+7 

"components”  with  values  of  +  10—  , 

The  compilation  of  the  program  requires  about 
4500  memory  locations,  excluding  the  array  storage. 


rv'*  •  ■ 

IMIS  ‘'It  Mtstr-^A',  aIlL  PH^'UilCf  A  C  i»L<()«-'‘.c*»sT  Sr>iTMt‘'l^ 

►  lALK'**  tT  Al.  ItFt  Man,  CAl  IHtOHY.  Vf>L  CT-11  NO  I  "ABlM  l4s4, 

►  ROM  A  GIVFN  DC  lUANbctR  PJNsTIOn  VFltOH  ITS  OHiVTMG  POlMT 

AOmiTANC»:.  THt  RESIST. *RS  urAIOWrO  ARE  NiNRtRtO  SERUfNTl  AtL  t  . 


IHOS  NOlt  0  IS  n.  node  l  is  l,  2  is  2»  li  is  5*  l?  is  If  21  IS  s, 

22  IS  bt  111  IS  NUMHtREn  T»  ETC.  THE  SUrtSCRIPI  0  MEFfHS  TO  GROIIN'T. 

% 


real  RET.SuM.MxA.MNR.xXR.niPP.BESTLAMY 


COMMENT 


% 


real  LAM.n£L*<'UMr»SUMC»OtNUM*NOMERl#NllMEM2A 

INTEGER  il*P»V,RtR.l»J»K»LtX»T»TOP.OOHF«  UEGREE*NJRlNP'lTS.MNP*MXP.WXKt 
FORMAT  rRMTn('OEGREC:**I.T*X5«'NO.  iNPl iTSS*  •  I i* A1 .0 ) « 

FORMAT  FRMTII  * ThE  DENOMINATOR  COrFFlrlENTS  ARE*.A2,0. 

ISR1<».S»A1  .Oils 

format  FRMT2(<THe  NUMERATOR  COEPFICIEmTS  OF  £nTRT**13«'  ARE**A2.0» 
(SR1«».A*A1  .fl)  )» 

FORMAT  FRMTSOTHE  NUMERATOR  COEFFICIENTS  OF  TI«*ol  ARr«»A2.0* 

<SR1M.S»AT .0) )» 

OM.T  FRMTRCTHE  denominator  COEFFICIENTS  OF  VIotOl  ARE**A2.0» 

IlM.BtAl .0) )» 

>  j  "aT  FRMT5(T3*X2»13tXR**N0Nf**Kl3»*M0NE*.Al.0)» 

►  ORMAT  FrMTA(I3»X2*I3»X3»R1<I.A»X«**NOmE**A1.0I» 

FORMAT  FrMT7(T3*X2«13*X6**N0nF**XR*R1U«R* Al.OIS 
format  FRMTR(T3tX2»I3*X3.RlR.P*X3*Rl4.a*Al.0lt 

FORMAT  FrMT9(«THE  FOLLOWING  COMPONENTS  ARE  iN  FArAQS  AND  0HMS«*A2.U» 

•THE  FRENI tENCY  IS  ONE  RAOTAN/SECOaiO*  » A1 .0* 

•THE  TREE  components  ARE^.AI.O* 
•NOOE**Xd.^CAP^*X13»^NOOe^*X8**RES*»A1.0l« 

FORMAT  FRMT1U(13»X3*R1R.'1»X8«T3*X3*R1II*R*A1.0I* 

FORMAT  FRMTIIMThE  TERMInATInO  COMPONFNTS  ARE**A2.0» 

•node  lNP|iT^*X6*^CAP^*XlM>*RES**A1.0)S 
FORMAT  FrMTUIEOI* 

FORMAT  FRMT 131 •VECTOR  NUMERATOR  COEFFICIENT  MATRIX^ * A».0. 

•NODE  INPi>T^*Al.UIS 

FORMAT  FrMTIMI  13»X3*13*XMf  (9ri2.S*A1  .OHA 
format  FRMT15(X6>I3*A«tl9Rl2.A«A1.0))V 

FORMAT  FRMT16I •VECTOR  OCNOMlNaTOR  COEFFICIENT  MATRIX^ .A2.n* 

XA* •nOOE**A1*0)A 

FORMAT  FRMTl 71 •DRIVING  PUInT  ADMITTANFC  DENOMINATOR  COEFFICIENT  MATRIX^* 
A2»n»X6*^A'0n£^*A|.0lft 

COMMENT  % 

sTRT..RCAO(0£FREE*NUMtNPUTSl» 

mRITEIFRmTI?)* 

R=DE6RECb 

4SNUMINPUTSS 

rOPSl»P*RU2» 

UONCsU**(R>in-2« 

mEGIN 

real  array  Ato..RtO..TOP.O..MIt  H»OtU. .RtO. . TOP) *CAP* OESIO. • TOPI > 


2Y.2C<OONF..rpP»U..Ml« 

mEAD  data  anu  print  it 


COMMrNT 

% 


READIFUR  Js(ltltM)  DO  FOm  AS(Otl*R)  DO  A(X>0»J}» 

FOR  JSIO.I.R)  DU  B<.i*0)> 

FOR  JslOil.R'll  00  tlJ»U))S 
MRirEIFR«Tl*Fr«R  k:(0*1»R>  jO  P(R*n)l» 

FOR  JxdtltMl  00  HRlTE(FR4r2« JaFOR  K  =  IU*1«R)  do  a(K»O.J))A 
«RITilFRMT3*F''R  KS(0*1*RI  oO  ■(K*n))» 

«R1TEIFRmT<i»F0R  XS(0*1»R-1I  uO  0(M*0))t 

CO'MMFNT 

sTARI  THc  synthesis 


« 


FOR  UX(0t)|«2*nuNcl  nu  BEv*lM  FOR  Pr<U/»*lttM  JO  BEGIN 
►  OR  ixiOtl.RI  JO  FOR  Ou  At  I *PtUl =A( I *Pt n|  44  1 1 .P. J) A 


cALfOLATt.  l»»Lf  rtt^T  1M| /1N>,  LA*l’<OA  AMD  CHOOtit  CLOSbcbT  ONE  TO  1/2 

% 

►  Oh  1)0  Eup  nO  rt^OlN 

kESTLAMSo  % 

oENOM:u(k«PI-n  (n«P) •(>(»•.< I /0(n»P)>e(H*P)«0(X-l  «P|/U<P-ltP)% 
nUMCM1SA(X*P*  Y)>HlU*P)*r)lX*P)/U(U*P>S 
NUMEhJSAI  X»o»  v»*o(0*p»  *01  X»P1  /[)IO»P)-p(X»P»» 

IF  otMOM  FOL  f  Then  i»n  To  jUmpoot  » 

t.AM=N(IMEHl/OcMUM  * 

.»EL=-N0MtP2/L)FN0M  i 

IF  (lam  toPu  1)  OW  (LA'*  LtO  0»  THEN  LAm  =0  % 

IF  (i»FL  »Fw  l»  OH  (DEL  LtO  0>  THEN  OEl  =U  t 
IF  AH<;(r)EL>n*'^)  LSS  AOi(LAM-O.b)  THEN  LAMsoEL  S 
IF  Ah^(La'4>0*<v)  LSS  AMS(hESTlAM-(J.5)  THFN  RFiTLAM=LAM  % 
jUmPUIM  .  . 

end  ptck  of  best  LAMhOA  > 
iF  resTlaM  eul  n  Then  bf>TlAm  =o.r  » 

COMMfNT 

HE  move  rc  components  ano  Cal  f  nek  aohittances 

% 

lAmshpSTlAm  * 

.»El=1-BESTlAH  s 
B(U»2«P«|l=R(''*Plt 

►OP  I=(l*l»p-1)  00  B(T*?*Pa1)=B(1*P)*LAM 

♦(»l(U«Pl*0(|»P)A0EL/0(o»p))  -  (q(R*P>*0(I*ltP)ALAN/o(R-l»P))» 
CAP(2«PA|):q(B*l*2*P«1)/o(P*l»P)ft 

FOP  l:(o»l*P-P>  1)0  0(I*?«o>lir0(l»PT-Mn«2AP«1)/CAP(2*P«i)  $ 

FOP  l  =  (n*l*P<-p)  00  R(I*?*PA2lrB(l»l*P}«n£L 

♦  B(H»P>a0(T*P)AlAm/0(h>1*P)  -  n(0*PlA0(  l*l*PUr)EL/''(0»P)ft 
h(P»1»2*pa2>=C‘(P»P)S 
HES(2«PA^)=*)(0>P)/B(Uf2A»>A2>S 

FOP  lr(n*l*P>P)  00  0(I*2*PF2)=0(lFl*P)>PES(2«P>2)AB<lAlr2APA2)  % 

Comment 

cAlCulaYe  nE«  tpansfep  Function  niimepatoh 

% 

FOP  lr(t*l*P*1)  OO  PEOlN 
MXA=U  » 

FOP  A  =  (1*1*M)  00  MXA=MAX(MXA*A(I*P*Kn« 

FOP  A:(1*1»M>  00  IF  A(1iH*a)  FBL  MXA  THEN  mXK=A  « 

MXBSMAX(B(lt2*PF|) tB(I*l*2«PT2) )« 

MNBSMTN(h( l*2APAl)*e(T*l*»APT2) )* 

IF  MXA  E«L  AIT-i»2*PT2)  IHEN  WXPSPPPT?  ELSE  MXP  s  2*P*l  » 

IF  MNP  EoL  R(T-l»?*PFP)  THEN  MNPsPAPT?  ELSE  MNP  =  2*Pf1  \ 

sUm:A(1»P*0)  « 

IF  S0'«  LtO  MNo  then  hFGIk 

POP  K=(1»1«M)  DU  BEGIN  A(l*MNP*A)SA(I.P*X)t  A(I>vXP»Kl:o  %  ENOS 
00  10  XFFHOONE  »  •‘NOS 

|F  SUM  LtO  MXO  Then  hFGI.) 

FOP  K=(1*1«M>  ou  BEGIN  A(I*MXP*R)SA(I»P*X)S  A ( I » uNP*X ) :J  S  ENQS 
-.O  TO  XFEPOUNE  i  ENOS 

FOP  j:(t*l*M)  00  «£6iN 

TF  •(I*P.JI  ruL  '1NA  Then  begin 

FOP  do  BF<*IN  A(l»MXP»R)=A(ItPft()S  AII>mnPiK)=U  S  ENOS 

A(I*MNP» J)=MNH  % 

AlI*M||P» J)=U  t 

ftU  TO  XFEPOONE  *  FNOS 

TF  A(I*P»JI  EGL  mxr  ThCn  begin 

rOP  KSlltliMl  OU  BFxiN  A ( 1 *MNP*A ) SA ( I »P»h )S  Al I » MXPt K ) sj  S  ENOS 
AdtMkP* J)=MXR  % 

A(I>MnP*J)=U  % 

r,0  10  XFFPnuNE  %  PNDS 

CnO* 

|F  M<«  GIP  MNP  then  hTgI  • 

rOP  A  =  (|«1*M)  ou  BE»T  4  m(1»MXP*RTSA(NPm)S  AI  I  •  ^nP^m  IFv)  S  CNQS 
a(I*mnp»vkx)  :  MNd  •> 

A  (  I  »  H  )  :«|XM»M'Im  % 

r.o  10  XPr^rUjNt*  CNOk 


»F  (sIM  v,T^  A.-JO  (•IK*  i.s^  T  IFN 

<0 

FU9  K  =  (  1  .  1  nu  UE  il>4 

IE  (<0M  ♦  AdfP.KlI  r.LO  MNR  THEN  <«EbIN 
T|FEr«IN.<  -  illM  % 

FOR  JS(1«1*F-1I  r)U  rtF^lnf  A<I*M*lP*j):A<t>Pf J|tA(I*MXP,j):ntE,4nt 
A(I»mnP»KI  :  OIFF  % 

Al  Ir*«|P*A):A(l»p»A)  .  OlFf  % 

FOR  00  riFblN  A(I*MXP»J»:A(l»P»J)«A(l«*«NP»J>=n$ENnS 

60  ro  XFEAOUMi  k  EMOA 

<;u'4  =  SUM  ♦  A(IfPtK)  » 

£Ni»  A 

e-'w  s 

XFCRO'VNF.  . 

ENO  transfer  function  CAlCUCATIOn  % 

FOR  A^dtltNI  00  FOR  i:d»l*Nl  00  A(  l.l*2ApA2*K)sA  ( I  •  9«Pk?*K  I  » 

FOR  k.sll«l*A)  00  AE61N  A(0*2RPAltK):A(U*P*A)S 
A(M«lf2*PA2*K)=A(R*P*K)«  CnOS 
tNO  P  » 
r:r-1  % 

i.HO  0,  RENOVA)  of  RC  TREE  % 

FlNACSTNIHESlR.. 

«RlTeiFRNT9)» 

FOR  PsIOdtOONE)  00  «RlTC(FllMnO»2APAt  »CAPI2*P«1  >  *2«Pf2*hES<2*Pa2M« 

CONmFNT 

FINISH  STNTNERIS  NOV  THAT  TRANSFER  FUNCTIONS  ARE  OF  UnITT  OEOPEE 

% 

riRlTtlFRATlDR 

FOR  US(OUNCA1*I*TOPI  OO  .^EOIh 

nURVpoS 

SUNCsnS 

FOR  00  SUMTsStMT  ♦  A(lki)»J)« 

FOR  J=(1*1*VI  00  SUNCsSimC  F  A(0»ii*J)« 

2T4tl»n):(R(lfli»-SUNV)/0(u*U>  % 

/C«u»n)s(mo*ii)>sui(C)/o(u*u)  « 

FOR  00  ZTlU*J>sAll*u*J>/0(0»UI  % 

FOR  Js<l«t*V>  00  ZC(U«J)sAie*(l*JI/0(U«U)  % 

FOR  j:(0tt*NI  00  BEOlN 

IF  2Y(u».t}  eoc  b  Then  seoin 

IF  7C(U*J)  rOL  0  ThfN  RFOIN 

NRITE«FRNT5*U*J|»  60  To  NEXT  % 

CNO  » 

vR|TC<FRHTb*U*j*2C<UtJ)IS  «0  TO  NEXT 

CNO* 

IF  /C(J»Jt  feOc  U  then  8E61N 

«R1TCIFRNT7*U*J»  1/FV<U»UM« 

60  TO  next  Enos 

«RirE4FHMTA*U*J»7C(.tf J)»  l/2Y(U*.in« 

<4Ert..  End  * 

cNO  transfer  rtnthesis  « 

•R1TC(FRnT12)« 

CONNRNT 

INSERT  PrInTEO  CIRCUIT  RuARO  PRINTOUT  SU8PR06RAN  HERE 

6 


CONHfnT 

PRINT  Out  COEFFICIENT  MATRlCt^ 

« 

RPOEbHE  % 
vRlTttFRvTlE)* 

«RlTt(FR4TlR)T 

FOR  u:(0*<l42*TOPI  00  RE61N 

FOR  P:IU/^*ltUt  00  «Et»Ii4 
A  =  l  » 

RRiTC<FRNTl<t*P«X*FuR  IStUti.H)  00  A(I«P*A)>« 

FOH  ik:(F*l*Rl  .tO  VHTTCtFRMTlSkAkFOR  1:(0»I»R)  uO  A(I*P.M|I$ 
F^OS 

RS9.I  % 

rp^n 


«RIUIFRttTl2)« 

FOR  U=(0*U^?f TOP)  00  RtGIM 

FOR  P:(U/2»i»U)  00  «R1TE (FRmT15*P*P0R  I=(Ofl»R>  OO  R(I*P))« 
P=R-l  % 

FND» 

RSOEGOCC  « 

•RITE(FRMT12)« 

mriteifrmtitu 

FOR  U=(Ofli«?*TOP)  00  RE61N 

FOR  PrfU/2»i*U)  00  if<RlTE(FRMT15tP»F0R  Is(0»l»R-l)  00  0<I»P))S 
R=R-1  % 

ENOS 

«RITe(FRNTl2)< 

END  PMOLE  thing  S 


COMMENT 

THIS  completes  the  SYNTHESIS  OF  ONE  VFCTOR  OF  THE  MATRIX. 
GO  BACK  AND  00  THE  NEXT  VECTOR. 

% 

GO  TO  SThT  S 
FINISH  t 
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THL  0CN0MIi^AT0»4  COf.FFICILNiS  AHF 
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As  an  experimental  verification  of  the  computer 
program  a  double  notch  filter  was  designed  using  the 
normalized  transfer  function 

u  -  _ <s^  *  1)  4) - 

(s^  +  0.04s  +  1) (s^  +  0.08s  +  4) 


which  produces  a  notch  at  j  and  at  2j, 

7 

Following  Hazony  and  Joseph  this  can  be  re¬ 
written  as 

U  =  ^  ^32 _ 

-  B(l-A  t3^) 

4  3  2 

Choosing  A  »  1  and  B  =  s  +  9s'^  23s  +  18s  +  4 

an  RC  transfer  function  vector  is  obtained  which  is 
almost  identical  to  that  of  Example  2.2; 

^  _  8.88s^-H7.9968s^+17.76s  ^  s'*->-5s^-*-4 _ 

”  s'*+9s^+23s^+18s+4  s^^s^+23s^+18s->-4 

The  driving  point  impedance  was  chosen  to  be 

V  s  s'*  9s^  -I-  23s^  +  18s  4 

^33  3  2 

4s  +  27s  +  46s  +  18 


This  network  was  synthesized  by  the  computer 
and  in  a  subprogram  the  notch  frequencies  were 

shifted  to  60  and  120  cps  and  the  component  values  < 

were  magnitude  scaled  by  a  factor  of  1.5  X  lO”^. 

i 

I 
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In  constructing  the  network  the  resistors  and 

capacitors  were  chosen  with  a  1  %  tolerance  and  the 

unity  gain  amplifier  used  was  a  D.C.  emitter  follower 

circuit  with  A  =  0,998,  Z.  =30  Meg,  and  Z  .  =  1  k. 

in  out 

The  response  of  this  circuit  is  shown  in  Figure  A. 2 — 
some  60  cps  noise  hampered  the  response  at  that 
frequency. 


Gain 

(db) 


I 

o 


Fig.  A. 
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2  Experimental  transfer  response  of  the 
double  notch  filter. 
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